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Abstract 

A method to precisely calibrate the oscillation amplitude in Dynamic Scanning Force Microscopy 
(DSFM) is described. It is experimentally shown that a typical electronics used to process the dynamic 
motion of the cantilever can be adjusted to transfer the thermal noise of the cantilever motion from 
its resonance frequency to a much lower frequency within the typical bandwidth of the corresponding 
electronics. Therefore, the thermal noise measured in the in-phase ("phase") and out-of-phase ("ampli- 
tude") output of such an electronics can be related to the thermal energy kT. If the force constant of 
the cantilever is known then the oscillation amplitude can be precisely calibrated from the thermal power 
measured in these signals. Based on this concept, two procedures for the calibration of the oscillation 
amplitude are proposed. One is based on simple calculation of the Root Mean Square (RMS) measured 
at the outputs of the electronics used to process the dynamic motion of the cantilever, and the other 
one is based on analysis of the corresponding spectrum and the calculation of the quality factor, the 
resonance frequency and the signal strength. 
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I. INTRODUCTION 



Scanning Force Microscopy (SFM) is an extremely versatile tool for Nanoscience that has a 
very high resolution, but usually only a modest precision. When SFM is used as a microscope, 
precise calibration of piezo movements is necessary, as well as careful control of piezo hysteresis 
and non-linearities. When used as a force measuring instrument a precise calibration of the 
displacement detector as well as the force constant of the cantilever is required. Dynamic SFM 
(DSFM) is by now probably the most extended SFM mode, since it allows operation at very 
low forces, and operates in the non-contact regime [1—4]. In DSFM the frequency (or phase) of 
the tip-sample system and the oscillation amplitude are the basic signals that are measured and 
should thus be calibrated for precise measurements. Precise measurement of frequency is easy and 
is implemented directly in most instruments. Calibration of oscillation amplitude is an issue that 
has received little attention, and less than what we believe it deserves. In Frequency Modulation 
DSFM (FM-DSFM) the frequency shift is related to (conservative) tip-sample interaction and the 
oscillation amplitude to dissipation. For Amplitude Modulation (AM-DSFM or "tapping mode" ) 
this separation is less evident, nevertheless the oscillation amplitude is an important parameter, 
since it critically determines non-linearity of tip-sample interaction [5, 6], onset of bi-stability 
and chaos [ )] as well as the precise interaction regime (attractive vs. repulsive) used for imaging 
[10]. Finally, as shown recently, if correctly calibrated the oscillation amplitude can be used to 
obtain a reliable image of the "true topography" that is independent of feedback parameters and 
scanning speed [11]. In the present work we will discuss a method to precisely calibrate the 
oscillation amplitude in DSFM. Essentially, the method relates the thermal noise measured in the 
"amplitude" and "phase" outputs of the electronics to the thermal energy kT. If the force constant 
of the cantilever is known then the oscillation amplitude can be calibrated from the thermal power 
measured in these signals. 

II. CALIBRATION PROCEDURES 

In the present work we discuss the response of the DSFM Dynamic Unit (DSFM-DU) in the 
frequency domain because we will analyze the response of the system to thermal noise Zth(t). Un- 
fortunately, thermal noise is non-trivial to describe using a "coherent" superposition of harmonic 
signals, as it is implicitly assumed in the Fourier Transform z(t) = J dv z(u)e l27Tiyt . The correct 
description of noise in the time domain should involve wavelet transforms [12], which is however 
beyond of the scope of the present work. We note that in the context of the present work the 
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DSFM-DU is used in a non-standard way, where no oscillation is applied to excite the cantilever, 
and the reference oscillator frequency v r is not at the resonant frequency of the cantilever. In- 
stead, the DSFM-DU is adjusted so that v r v — bw/2 where bw = 1/r (r: time constant) is 
the bandwidth of the DSFM-DU. As discussed in detail in Appendix A. The thermal peak of the 
cantilever appears centered within the spectral range of the output of the DSFM-DU. 

Essentially, the reason for this frequency shift is that when the normal force signal is internally 
multiplied with the reference frequency v r of the DSFM-DU, frequency components at Vq — v r 
and vq + v r appear (for example from terms like cos(27rz/ t) cos(27rz/ r t), see also figure 1 for a 
graphical representation). Usually, the DSFM-DU is opperated with v = v r , then thermal noise 
is observed at very low frequencies v ~ 0, that is, at DC. At this low frequency the thermal noise 
is usually not clearly distinguised from pink noise (1/V-noise) of typical electronic components. 
However, when vq ^ v r , the thermal peak is shifted to the difference frequency vq — v r , which 
may be well below the resonance frequency, v . As discussed above, we typically adjust v r so that 
Vq — v r w bw/2 « 4 kHz in our system. 

Figure 2 shows experimental data for the deflection signal Ud(t), and the output u y (t) measured 
when the DSFM-DU is configured as just discussed. The main graphs show spectra, and the insets 
time domain data. The spectra have been acquired using a Spectrum Analyzer [13] as well as 
the "noise" function of a commercial Lock-In Amplifier [ I I] . The lines through the data points 
correspond to a "Lorentzian + offset" fit function f{y) [15, 16], 

/M = (i - w»Q)r +el= 4Jg(l/)|2 + e - (1) 

where g(v) is the mechanical gain of the cantilever (see relation (Al)). The corresponding pa- 
rameters for the different fits are specified in table I, and discussed in more detail below. The 
scattered points around the horizontal lines show the error of the experimental power spectrum to 
the corresponding fit, which show little tendency for both spectra, therefore we conclude that the 
chosen fit functions describe well the experimental data. Both spectra lead to essentially the same 
results, in particular for the values of the resonance frequency v Q and the quality factor Q. Since 
the DSFM-DU internally amplifies the normal force signal, the signal strength of the deflection 
signal Ud{t) is smaller than that of the DSFM-DU outputs u x (t) and u y (t), leading to different 
thermal noise e t h and electronic noise e n for the deflection spectrum Udiy) compared to u x (u) 
and Uyiv). We therefore conclude that the outputs u x (t) and u y (t) of the DSFM-DU essentially 
reproduce the thermal noise spectrum of the deflection signal and that the DSFM-DU therefore 
"sees" (internally amplified) thermal noise shifted to the frequency v Q — v r . It should be therefore 
possible to apply the well-known thermal calibration method based on the equipartition theorem 



[16-22], 

|<4(*)U = \kT (2) 

to the calibration of the oscillation amplitude. 

We note that, as for the case of normal force calibration, relation (2) has implicitly two un- 
knowns: the force constant c of the cantilever and a sensitivity calibration /3 (units: nm/V) to 
convert the measured signal u(t) (Volts) into physical length units (Nanometers): 

z(t) = 0u(t) (3) 

In the present work, we will assume that the force constant is known or has been determined 
independently (see, for example [22, 23]), and will thus use (2) to calibrate the sensitivity (3 of 
the DSFM-DU. We will present and discuss two procedures based on equation (2): a simple one 
which only requires estimation of the RMS noise of the output of the DSFM-DU, and a second 
one taking into account the full spectral response of the outputs u x (t) and u y (t). 



A. Calibration based on root mean square estimation 

The first calibration procedure is implemented experimentally by acquiring the outputs u x (t) 
and u y (t) for a certain period of time with the laser on and with the laser off. The noise power 
that is related to thermal noise is then, 

{u th f = (u 2 on (t)) rms - (u 2 oS (t)) Ima (4) 

The reason to acquire data with the laser off is to estimate and subtract electronic noise that is 
not related to the thermal noise of the cantilever motion. In particular for hard cantilevers, where 
thermal noise is quite small (see below) the calibration may be severely wrong if this electronic 
noise is not subtracted. Without substraction, it is implicitly assumed that all measured noise 
(that is, also electronic noise) is thermal noise. The voltage noise Au^ h = J ((u^ h ) 2 ) of the 
output u x (t) is related to the "phase" noise while the noise Ait^ = J {{u v th f) is related to the 
oscillation amplitude. Relation (4) can of course also be applied to the fluctuations measured 
in the deflection signal Ud{t), then Auf h = J ((u^) 2 ) is the deflection noise (in Volts) of the 
cantilever motion. We recall that thermal noise is distributed equally between the in-phase and 
out-of-phase components; then, if the amplification of the outputs u x (t) and u y (t) as well as the 
deflection Ud(t) are the same (the latter is usually not the case, see below), 

Au th = A< = Au d th /V2 (5) 
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since the total deflection signal and the in-phase and out-out-phase componentes x and y are 
related by u 2 d {t) = u 2 x {t) + u 2 (t) and therefore f3 x = (3 y (see, for example, [15]). The values 
Awy, and Au y h can be either computed directly from the time domain signals u x {t) and u y (t), 
or by adding (integrating) the noise spectrum in the frequency domain. Due to the different 
normalizations found for the Fourier Transform and the Power Spectral Density in the literature 
as well as in software packets, it is usually more secure to use the first option. The insets presented 
in figure 2 show the data used for this calibration procedure. The data from the DSFM-DU is 
presented as scatter plot (u x (t),u y (t)); as it would be visualized with an oscilloscope in x — y 
mode. Data is shown for the laser on and off. The diameter of this scatter plot is proportional 
to the fluctuations in both directions. The data corresponding to the deflection signal Ud(f) is 
shown "normally" as a function of time, again data is shown for the laser on and off. Calibration 
factors are computed directly as (3 y = z t h/ Au v th / \/2 (see relation (5) for the factor y/2) and 
j3 d = z t h/ Auf h with z t h = \JkT/c rms thermal noise movement (unit: nm, from equation (2)). 
The results obtained with this first method are summarized in the left column of table I (termed 
"first procedure"). The calibration factors for the cantilever deflection and for the (amplitude) 
output of the DSFM-DU have been measured for a cantilever with c = 1.6 N/m, as determined 
using Sader's method [23] from the resonance frequency and quality factor (see table I). 

B. Calibration using the thermal noise spectrum 

The second calibration procedure is based on estimating the area under the thermal noise 
spectrum using the values determined from the fit to a "Lorentzian function + electronic noise 
function", f{y) from relation (1). In this context we recall that J °° dv |g(^)| 2 = Q Vq 7r/2, 
therefore the measured "thermal noise power" as determined from the fit parameters is e 2 h Q uqti/2, 
where e 2 h is the strength of the thermal signal (units: V 2 /Hz). As discussed previously, in the 
frequency domain the DSFM-DU shifts the original deflection signal z(u) from a signal centered 
around the resonance frequency z/ to signals u x {y) and u y {y) centered at u — v r « bw/2 where 
v r is the reference frequency adjusted by the DSFM-DU. In addition, the original signal Udit) 
is amplified by the internal gains of the DSFM-DU. In principle, there are two possible ways of 
processing the data obtained from the DSFM-DU: the spectra |u x (z/)| 2 and \u y {v)\ 2 may be fitted 
either directly to the function f{y) from relation (1), resulting in a resonance peak near bw/2 
(see figure 2) -or the frequency shift induced by the DSFM-DU may be "undone" by adding the 
reference frequency to the frequency of the power spectra u x {y) and u y {y) in order to fit these 
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back-shifted power spectra u x {y + v r ) and u y {y + z/ r ), 



• Option 1: u x (t) raw data — > FT and PSD — > \u x (u)\ 2 — > fit — > {e t h, Q, vo,e 
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} 



• Option 2: raw data — >FT and PSD— > |-u.j.(z/)| 2 — > freq. shift — > \u x (u + u r )\ 2 — > fit 



->■ {e t ft, Q, Vq, e n } 

The second option is the correct one, as relation (2) holds for the "true" deflection signal with the 
thermal noise peak at the "correct" resonance frequency vq. 

We therefore fit the frequency shifted spectrum \u x (v + v r )\ 2 to the "Lorentzian + offset" 
function f{v) from relation (1), in order to obtain the "strength" eth of the Lorentz function, the 
quality factor Q, the resonance frequency z/o and the noise power e n that is not Lorentzian. For 
the experimental spectra of the output Ud{t), relation (2) is rewritten as 



where (3 d (unit: nm/ Volts) is the factor that converts the deflection signal Ud(t) (in Volts) to the 
physical amplitude z{t) (in nm): z{t) = f5 d Ud{t). Similarly, f3 x and (3 y convert the signals u x (t) 
and u y {t) into the in-phase and out-of-phase components of cantilever deflection. For the outputs 
u x (t) (and u y {t)) relation (2) is 



The last relation for f3 y in (7) is obtained because in our case the amplification of the two outputs 
u x (t) and u y (t) is the same; therefore, as discussed above, (3 X = (3 y . The parameters Q, u and 
e t h are obtained from the fits, kT is the temperature of the cantilever, and c = 1.6 N/m is the 
force constant, determined by Sader's method [23] from the Q factor and the resonance frequency. 
The right columns of table I (termed "second procedure") summarize the results obtained from 
this second calibration procedure when applied to the deflection data Ud{t) and to the out-of- 
phase output u y (t) of the DSFM-DU. We will now comment the different fields of this table. 
The Lock-In amplifier is used to measure the noise density, and does not acquire time domain 
data, therefore the first procedure cannot be applied, and the corresponding fields are empty. 
The Spectrum Analyzer used allows to simultaneously measure the real-time data, as well as the 
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spectrum of this data, therefore the first and second procedure can be applied. For this cantilever 
with a relatively large thermal noise (about 50 pm) the first and second calibration procedure 
give similar results, although the second procedure seems to have less error (see discussion below). 
In order to compare the results obtained from the deflection signal with that of the output (s) of 
the DSFM-DU a normalized calibration factor [5^ has been introduced that takes into account 
the internal gain of the DSFM-DU. This normalized calibration factor should be the same for the 
deflection data and the outputs of the DSFM-DU, which is indeed the case within the experimental 
error of the measurements. 

Figure 3 shows the spectra of the deflection sensor and the outputs of the DSFM-DU for a 
hard cantilever (c = 67 N/m as determined by Sader's method [23]), with a quite low thermal 
noise amplitude (z t h = y/kT/c ~ 8 pm). In this case, data is shown for different gains of the 
DSFM-DU, and the spectra have been acquired directly from the (calculated) power spectrum 
of the time domain signals acquired by the SFM-control unit (more precisely: from ^(t) -images 
acquired with the tip far from the surface at maximum adquisition speed, without scanning and 
no excitation applied to the driving piezo of the cantilever). Note that for hard cantilevers, the 
resonance frequency may be easily outside the acquisition bandwidth of the analog to digital 
converters of the SFM-acquisition electronics (usually 16 bits or more and thus rather slow) and 
therefore the corresponding thermal noise of cantilever motion cannot be "directly seen" in the 
digitalized normal force (deflection) signal Ud(t). By contrast, the thermal noise peak is easily 
brought into the bandwidth of the analog to digital converters when the outputs u x (t) and u y (t) 
are used, since this noise is now at the much lower difference frequency. For this hard cantilever, 
the very small thermal noise signal is usually significantly smaller than the electronic noise or 
other non-thermal fluctuations. In our setup, this is indeed the case for the direct deflection 
signal as well as for the outputs of the DSFM-DU when its internal gain is smaller than 10 (see 
the "scatterplots" in the inset of figure 3). Then, it is essential to "normalize" the measured signal 
by subtracting the electronic noise (laser off) from the total noise (laser on) as discussed above 
(see relation (4)). Still, for low gains this procedure does not give satisfactory results; then only 
the second procedure method is precise. Surprisingly the second procedure still works for gains as 
small as g = 1 and for the direct deflection signal. Note that for these low gains almost all noise 
is electronic noise. The fit to the function f{y) = e^ h \g (z/)| 2 + is thus a very effective way of 
"filtering" all non-thermal noise. For high gains the calibration factors rise by an amount that is 
not compatible with the error of our measurements for both procedures. We believe that this is 
not a problem of the calibration procedures, instead we think that this is due to low-pass filtering 
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of our DSFM-DU that reduces the (nominal) signal strength and leads to a higher calibration 
value at high frequency (more nanometers of deflection are "needed" for 1 V signal). 



III. CONCLUSION 



We have presented two methods for calibration of the oscillation amplitude in DFSM , one based 
on simple calculation of the RMS value of the output of the DSFM-DU u y (t), and the other one 
based on analysis of the corresponding spectrum u y [y) and calculation of the parameters Q, u , 
eth and e n . From the results summarized in tables I and II we conclude that both methods give 
consistent results, even though the second method is considerably more precise and robust, in 
particular for the case of hard cantilevers, where thermal noise has a much smaller amplitude and 
the signal to noise ratio of thermal noise (which in this case is "good" signal) vs. other noise sources 
is much lower. The second method has several advantages: first, it explicitly "filters" thermal 
noise from other noise sources, since it will only take into account signal that has a Lorentzian 
shape, signal not compatible with this shape is taken into account by the constant factor e n and 
the corresponding signal power is rejected for the calculation of the amplitude sensitivity. Also, 
fitting of a spectrum with many data points in order to obtain the four parameters Q, vq, eth and 
e n results in effective data averaging, and thus additional improvement of the estimation of the 
thermal noise power. And finally, from a strictly theoretical point of view, the first method is not 
correct because it is only an approximation valid for high Q factors. In fact, the first method only 
"sees" thermal noise in the small bandwidth 1/r C t'o around the resonance peak, but not all the 
noise under the Lorentz function, and in particular not the noise in the low frequency "tail" (from 
DC to i>Q — bw/2). This noise is "filtered away" by the DSFM-DU and thus not taken into account; 
therefore the first method underestimates noise, and overestimates the amplitude sensitivity (see 
equation (2)). Since the relation of thermal noise in the resonance peak to that in the low frequency 
"tail" is Q : 1 (see for example [15], section 3) this error is negligible for experiments in air and 
vacuum, but is expected to be significant for the low Q factors encountered in liquids. Finally, 
we note that determination of the parameters eth, Q and vq using the outputs of the DSFM-DU 
involves significant improvement of signal because of the internal gains of the DSFM-DU and the 
principle of Lock-in detection, which implies an important reduction of bandwidth. 
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Appendix A: Processing of thermal noise by a Dynamic Unit 

We recall that a typical electronics used for analyzing the dynamics of the cantilever (see figure 
1), may be implemented using a Lock-in detection scheme [1, 2, 15]: the input signal Ud(t) to 
be analyzed is multiplied by two reference signals in quadrature (a r cos(2ni' r t) and a r sin(27rz/ r t)) 
and then filtered with an appropriate time constant r. The corresponding output signals of the 
DSFM-DU are the in-phase ("phase") and the out-of-phase ("amplitude") signals 1 u x (t) and 
u y (t). Generally in DSFM the reference signal a r cos(27rz/ r t) is used to excite the cantilever at 
resonance, then u x = 0; and changes of u x (t) are proportional to the phase (and thus frequency) 
variation of the cantilever oscillation and u y (t) is its oscillation amplitude. For FM-DSFM a 
feedback loop (Phase Locked Loop, PLL) is used to track the resonance frequency of the tip- 
sample system by changing the excitation (i.e., the reference) frequency, in order to keep u x (t) = 
(not shown in figure 1). In a recent work [15] we have analyzed in detail how such a DSFM-DU 
processes signals in the presence of thermal noise. In particular, it was shown that the DSFM-DU 
"frequency-shifts" a signal Ua(t) = Re[a(z/)e 27rwt ] with 

a{y) = x{v) + iy{v) = - 2 ao = a g (v) (Al) 

1 - {u/U Q ) +1{V/V )/Q 

at its input, to the sum and difference frequencies, where g(y) is the complex "mechanical gain" of 
the cantilever. The outputs "phase" and "amplitude" have frequency components at z/ s = u + v r 
and i>a = v o — v r : u x {t) = Re [M(i)] and u y (t) = Im [M(t)] with 

1 (A2) 



M(t) = 



_i e -2iri(t vs+1/2) _^ \ c +2Tri(t wa+1/2) 



1 + il-KV^T 1 + iI'KVaT 

1 Note that in [15] the nomenclature (x q (t)) T and (y q (t)) T was used, that is: u x (t) = (x q (t)) T and u y (t) = (y q (t)) T . 
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where, the matrix notation in [15], appendix A, has been translated into the more compact complex 
notation. When the input signal is not at a well defined frequency, but distributed around a central 
frequency u , Udit) = f du Udiy)e l2lvvt , then the outputs of the DSFM-DU will be: 



where we have assumed that the term with the sum frequency z/^ (see relation (A2)) can be 
neglected because the time constant r of the DSFM-DU is much larger than the time l/z/£ (see 
also figure 1). The DSFM-DU therefore shifts the spectrum of its input signal Ud{t) centered at u 
to a frequency z/a = z/o — v r of the output signals u x (t) and u y {t). That is, if the power spectrum 
of Ud(t) has a maximum at u , the output signals u x (t) and u y (t) will have maxima at the much 
lower frequency z/a = — Vr, as stated in the main text of the present work. 



u x (t) = -Re [D (t)) and u y (t) = -Im [D (t)} 



(A3) 



with 




(A4) 
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Appendix B: Figures and Tables 




Figure 1. Schematic description of a typical Lock-in type DSFM detection unit. The signal to be analyzed 
by the DSFM detection is assumed to be centered around some frequency vq. It enters the detection 
unit at the input "in" , is amplified by some factor g that can be generally selected and usually high-pass 
filtered (for simplicity the corresponding components are not shown) before it is multiplied with two 
reference signals in quadrature at some frequency f re f- After this multiplication, the signal is shifted to 
the frequencies vq — u re f and vq + f re f- The resulting signals are then low-pass filtered to remove the 
higher frequency component {i>q + ^rci)-, an d brought to the outputs u x (t) and u y (t). 
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Figure 2. Power spectra u y {y) (left curve) and u<i{v) (right curve; this curve is amplified to match the 
vertical range of the amplitude data) of the outputs u y (t) of the DFSM-DU unit and the deflection signal 
Ud(t). Data has been acquired with an audio analyzer (green, small points), directly by the AD-converters 
of our control electronics (red, larger points) and with a Lock-In amplifier (blue, the largest points). Left 
insets: phase vs. amplitude, (u x (t),Uy(t)) data acquired with an audio analyzer. Right inset: deflection 
vs. time signal, u^t). For both insets: dark blue, signals with the laser on; light blue, signals with the 
laser off. 
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4 316 v (kHz) 

Figure 3. Thermal noise data acquired for a hard cantilever (c = 67 N/m). Left side: power spectra 
of the outputs u y (t) of the DFSM-DU unit acquired with different input gains, (a) Power spectra for 
gains g = 1, 2, 3, 10, 30 and 100. Note that in this figure a logarithmic representation of the spectra has 
been choosen in order to compress the vertical scale to fit all spectra with very different gains. Lower 
right graph: Spectrum of the deflection vs. time signal Ud(t), measured with a digital oscilloscope. Inset: 
phase vs. amplitude scatter plots (u x (t),u y (t)) from the data that was used to compute the spectra u y {v) 
shown left; dark blue: signals with the laser on, light blue: signals with the laser off. The upper three 
scatter plots of the graph correspond to (u x (t),u y (t)) data for gains g =100, 30 and 10; the lower three 
scatter plots to data for gains g = 3, 2 and 1. 
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Table I. Results obtained with the two calibration procedures applied to the measured deflection and 
amplitude noise acquired from a cantilever with a force constant c = 1.6 N/m. 
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0.36 


0.30 


0.21 


40 


40 


316050(±59ppm) 


570(±13%) 


9.4(±9%) 


50(±14%) 


51(±14%) 


(Uy(t)) G = 10 


0.73 


0.38 


0.63 


13 


44 


316054(±6ppm) 


570(±1.3%) 


33(±1%) 


14.1(±1.7%) 


50.1(±1.7%) 


(uy(t)) G=30 


1.53 


0.60 


1.41 


5.6 


48 


316053(±0.4ppm) 


576(±0.67%) 


74.8(±0.56%) 


6.21(±0.9%) 


53.2(±0.9%) 


( u i/(*)}g=ioo 


6.24 


2.02 


5.91 


1.33 


53 


316055(±0.3ppm) 


563(±0.6%) 


316.2(±0.51%) 


1.49(±0.8%) 


52.9(±0.8%) 



Table II. Results obtained with the two calibration procedures for a hard cantilever (force constant c = 67 
N/m). 
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